Single-bit Feedback and Quantum Dynamical Decoupling 



Francesco Ticozzi ' Tj and Lorcnza Viola 2 -0 

1 Dipartimento di Ingegneria dell'Informazione, Universitd di Padova, via Gradenigo 6/B, 35131 Padova, Italy 
'Department of Physics and Astronomy, Dartmouth College, 6127 Wilder Laboratory, Hanover, NH 03755, USA 

(Dated: February 1, 2008) 

Synthesizing an effective identity evolution in a target system subjected to unwanted unitary or 
non-unitary dynamics is a fundamental task for both quantum control and quantum information 
processing applications. Here, we investigate how single-bit, discrete-time feedback capabilities 
may be exploited to enact or to enhance quantum procedures for effectively suppressing unwanted 
dynamics in a finite-dimensional open quantum system. An explicit characterization of the joint 
unitary propagators correctable by a single-bit feedback strategy for arbitrary evolution time is 
obtained. For a two-dimensional target system, we show how by appropriately combining quantum 
feedback with dynamical decoupling methods, concatenated feedback-decoupling schemes may be 
built, which can operate under relaxed control assumptions and can outperform purely closed-loop 
and open-loop protocols. 

PACS numbers: 03.67.Pp, 03.65.Yz, 05.40.Ca, 89.70.H-c 



I. INTRODUCTION 

The role of real-time measurement and information 
feedback in designing control strategies is well recognized 
for a wide range of applications in classical control theory 
Feedback in the quantum domain, however, unavoid- 
ably involves additional challenges due to the intrinsically 
probabilistic 0, 13 and, typically, destructive character 
of the quantum measurement process Q - leading to a 
host of delicate issues and added difficulties for many 
tasks of interest to the emerging field of quantum con- 
trol. As a result, real-time quantum feedback is nowadays 
actively investigated at both the theoretical and exper- 
imental level. In particular, continuous-time feedback 
strategies, based on indirect weak measurements, have 
been originally proposed for cavity QED systems [j| |(| , 
and have recently begun to be experimentally accessi- 
ble From a control-theoretic standpoint, a general 
analysis of strategies based on average measurements on 
quantum ensembles and feedback has been considered in 
[3 • Suggestive applications of continuous measurements 
and feedback in the context of quantum error correction 
have also been investigated [Toj . 

Our focus in this work is to begin exploring possible 
uses and advantages of discrete-time quantum feedback, 
based on strong projective measurements, in connection 
with a fundamental quantum control task: engineering a 
desired unitary operation on a target system, which may 
be generally interacting with an uncontrollable quantum 
environment. Thus, given a set of available control oper- 
ations, the goal is to design a protocol capable to enforce 
a desired net evolution on the system at a given final 
time, irrespective of the system's initial state. This prob- 
lem, unlike the simpler case of state preparation, forces 
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in general the introduction of indirect measurements on a 
suitable ancillary system. The first step in this direction 
is to employ the minimum physical amount of informa- 
tion we can gather, a single bit. In particular, in what 
follows, we aim to construct and test protocols for en- 
gineering the simplest unitary operation (the identity or 
no-op gate |llj ) by using, along or in conjunction with 
additional control schemes, a feedback loop built from 
basic, arbitrarily fast and accurate control operations. 

The analysis and techniques presented here exploit and 
extend ideas close to those introduced in a number of ear- 
lier works. In particular, our ancillary qubit system may 
be seen as a concrete instance of a quantum controller 
or ; more generally, a universal quantum interface 
Il4| . Control schemes exploiting auxiliary quantum 
probe systems and measurement capabilities have also 
been recently discussed in However, while the 

main emphasis of the above-mentioned literature is on es- 
tablishing general existential results on quantum accessi- 
bility and controllability, our interest here is on providing 
constructive, synthesis results instead. 

Two main families of protocols emerge: on one hand, 
feedback- enacted decoupling protocols - which, in the in- 
finitesimal time limit, directly relate to feedback imple- 
mentations of quantum dynamical symmetrization proce- 
dures |l7lll8| and, when appropriate conditions are met, 
allow to relax the fast control limit of open- loop schemes: 
on the other hand, feedback- concatenated decoupling pro- 
tocols - whereby the output of a primary open-loop dy- 
namical decoupling (DD) protocol serves as input to an 
outer feedback block, or viceversa. Interestingly, the 
analysis of the feedback-enacted procedures may be also 
regarded as stemming from the communication-oriented 
approach of |l9U20l]. or as a particular instance of a quan- 
tum error-correction procedure |2lj |. In a similar spirit, 
an interesting discussion on the ensuing trade-off between 
information gain and disturbance for different measure- 
ment strategies in a single-qubit setting can be found in 
|22j . From a technical point of view, our main results fur- 
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ther demonstrate the usefulness of linear algebraic meth- 
ods for finite-dimensional quantum control settings, along 
the lines introduced in |2Jj |. 

The content of the paper is organized as follows. Fol- 
lowing a presentation of the general control-theoretic set- 
ting in Section [HI we develop our main theoretical anal- 
ysis and results for a generic finite-dimensional target 
system in Section IIIII By restricting to the simple yet 
paradigmatic case of arbitrary state stabilization for a 
single qubit, in Section ll VI we describe various concate- 
nated feedback-decoupling protocols built, in particular, 
from both maximal and selective DD schemes 0, |^| ■ 
The performance of different schemes in the non-ideal 
limit of finite control time scales is numerically inves- 
tigated and compared in Section A discussion and 
concluding remarks follow. 



In the reminder of this section, we recall some basic 
results related to DD and symmctrization, referring the 
reader to 0, 0, 0, |2|| for more detailed discussions. 

DD strategies for open quantum systems are based on 
the idea of coherently modify the target dynamics by in- 
terspersing the natural evolution with (ideally) instanta- 
neous control manipulations acting on the system alone. 
Let H c (t) <g> Is denote the applied, semiclassical control 
Hamiltonian. DD is most conveniently described by ef- 
fecting the time-dependent change of basis in Hs induced 
by the control propagator U c (t) = Texp{— i J Q H c (s)ds} 
(transforming to the so-called "logical frame", and set- 
ting H = 1 henceforth). For deterministic DD, control 
sequences are additionally assumed to be periodic, mean- 
ing that U c (t + T c ) = U c (t) for a cycle time T c > 0. Then, 
the total transformed Hamiltonian reads: 



II. SYSTEM AND CONTROL SETTINGS 

We will consider a finite-dimensional quantum system 
S, coupled in general to an uncontrollable, possibly un- 
known quantum environment £ . Let 7is and Hs de- 
note the system and the environment Hilbert spaces, with 
dim(7i_4)= ds, dim(Ti.s)= dg, respectively. Assume that 
the joint system plus environment dynamics is driven by 
a time-independent Hamiltonian of the form: 



H(t) =^Ut{t)SiU c {t) ® Bi 



(2) 



i=0 



The associated logical propagator U(t), which coincides 
with the full Schrodingcr propagator at every instant 
Tn = NT C , TV £ N, may be conveniently expanded in 
the Magnus series |2q . 



U(T N ) 



— iHT., 



N 



(3) 



H tot = H s ®Is + Is®H s + HssY^ S i® B ^ (!) 



where H denotes the average Hamiltonian [2J|. Explic- 
itly, the lowest-order terms are given by 



where A: is a finite integer, H$, Hg, Hss are the system, 
environment, and interaction Hamiltonians, and Si, Hi 
are arbitrary Hcrmitian operators on Us and Tig, re- 
spectively. Without loss of generality, we may identify 
Sq = Is and take the Si to be traceless for i = 1, . . . , k. 

We will assume access to the following control re- 
sources and capabilities: 

• An additional two-level system A with Hilbert 
space Ha (an ancilla system, or a "quantum infor- 
mation channel" henceforth), able to be prepared 
in a given pure state |</>). We will assume A to have 
a decoherence time much longer than the relevant 
control time scale, so that it can be treated as a 
closed system, evolving under the trivial Hamilto- 
nian Ha = under ideal conditions. 

• Arbitrarily fast and strong unitary transformations 
on the system and the ancilla; 

• Fast conditional gates, i.e. unitary operations of the 
form C# s = |0)(0U®7s + |l)(lU®t^onW^®Ws 
- {|0), |1)} denoting an orthonormal basis in Ha- 

• Projective "Yes-No" measurements on the ancilla, 
corresponding to the projections on the reference 
basis subspaces. Suppose, moreover, that we can 
read, record, and use the classical information we 
obtain to influence the subsequent control strategy. 



H (0) 



1 

tT 



dsH(s) , 



(4) 



= - — 

2T r 



dsi / ds 2 [H{ S2 ),H{ Sl )}, (5) 
Jo 



c Jo 

with higher-order terms expressible as a power series in 
the control parameter T c . Given a finite evolution time 
T > 0, DD is implemented by considering a time-slicing 
into TV control cycles over [0, T], with T = T N = NT C . In 
the ideal limit of arbitrarily fast control, N — > oo, T c — > 0, 
it can be shown the higher order terms in (|SJ) are negli- 
gible in a suitable norm, and the dominant contribution 
to the logical propagator takes the simple form 



U(T N ) 



In practice, higher-order terms turn out to be critical 
in affecting the control performance when the ideal con- 
dition TV — > oo is relaxed. The ideal limit still proves, 
however, a useful starting point. Suppose we choose a 
piecewise-constant control propagator, 



U c (t) = Gj, jAt<t<(j + l)At, 



with j = 0, 



1, 1 < n g G IN, and A< = T c /n g . 



From Eq. Q , this control prescription gives the following 
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overall average Hamiltonian: 




By ensuring that the cyclicity requirement is fulfilled (im- 
plying that Ylj Gj — Is), the above Hamiltonian specifies 
the lowest-order stroboscopic evolution in the physical 
frame also. In practice, for a large class of DD schemes 
the control operations Gj form a discrete group Q, with 
n g = \Q\. In this case, it is possible to show that the 
quantum operation in parenthesis in Eq. © effectively 
projects onto system Hamiltonians which are invariant 
under Q - hence commute with the Gj themselves. 

In what follows, we will exploit DD to synthesize a 
vanishing average Hamiltonian, i.e. to "freeze" the dy- 
namics to first- (or higher) order in time. Let us recall 
two relevant DD schemes for a single two-level system, 
which will be useful in Section Hvl 

(i) First, suppose the control objective is to suppress 
evolution due to a known drift and/or a known cou- 
pling Hamiltonian which, without loss of generality, may 
be assumed to be along z. Thus, the non-trivial Si in 
(JTJ are proportional to the Pauli matrix a z . Then, DD 
over [0,T C ] is achieved by letting S evolve till T c /2, then 
applying an instantaneous transformation a such that 
aa z cr — —<7 Z (e.g., a = a x ), and applying a second 
identical transformation after another T c /2 to complete 
the cycle. The desired averaging is ensured by the prop- 
erty o~ z + ercr 2 crt = 0. To improve performance for finite 
T c , a symmetrized version of the above sequence (so- 
called Carr-Purcell DD) ensures that odd corrections in 
the Magnus series are automatically eliminated. This is 
attained with a slight re-arrangement of the control time 
intervals: Let iS evolve freely till t = T c /4, set the control 
propagator to a and let it constant until t = 3T C /A; then 
reset the control propagator to Is and let iS evolve for 
the remainder of the cycle. 

(ii) If no information is available about the drift and/or 
coupling Hamiltonian to be suppressed, a maximal DD 
strategy need to be employed. As shown in ^tJj cycling 
the control propagator through the elements of an or- 
thonormal basis for su(2), i.e. I , o~ x , o~ y , cr z , ensures that 
the resulting average Hamiltonian is always zero, reflect- 
ing the fact that the action of the Pauli group is irre- 
ducible. Interestingly, while the order in which the con- 
trol propagators are switched along the cycle is irrelevant 



to lowest-order averaging, different "control paths" over 
the same set may lead to different control performance for 
finite T c , due to their different influence on higher-order 
corrections. 

III. FEEDBACK-ENACTED DYNAMICAL 
DECOUPLING 

Our first goal is to analyze the possible advantages 
of using a classical, single-bit feedback strategy to en- 
gineer the identity evolution on S. It is apparent that 
projective measurements on S itself cannot be of any 
use toward realizing a unitary operation, as in general 
the pre-measurement state is irreversibly lost after the 
information is gathered. Thus, as anticipated, we will 
consider a combination of idealized, instantaneous uni- 
tary control actions on the pair S plus A, and projective 
measurements on A, to dynamically eliminate unwanted 
components. 

Assume that A is prepared in an initial state \<fi) = 
a\0) + b\l), with coefficients a, b to be determined. With- 
out loss of generality, we may also assume that the £ is 
initially in a pure, though unknown in general, state - 
denoted by Finally, let denote an arbitrary pure 
state of iS. The basic steps of a Feedback- Enacted DD 
protocol (FDD) may be described as follows: 

(I) Rapidly entangle S and A, by performing a condi- 
tional gate Us'- 

c$ s \4>M = a\om + b\i)u s m. (7) 

(II) Evolve freely, in the presence of f, up to a finite 
time T. The combined final state may be described as 
follows: 

UtoVse {a\Q)mO+h\l)U s \m)) = 
= a\0)Use\m) + b\l)U se (W» 10 • (8) 

(III) Proceed by undoing the preparation step, that is, 
effect the following conditional transformation: 

C\ ® I s (a\0)Uss\m) +b\l)U S£ (U s \^)M)) = 

a\0)Us £ \ip)\0 + b\l)(ul ® I £ )Uss(U s ® IeW\0- ( 9 ) 

The above equation makes it clear that, in order to ob- 
tain a balanced averaging of the two unitary components, 
equal weights on the two ancilla states arc needed. We 
thus set a = b = 1/ y/2 henceforth. 

(IV) Apply a Hadamard transformation on A, 
obtaining the combined state 



V2 



IS£ 



\0)UssMO + \l)(Us®l£)Uss(Us®Is)\m) 



\\0) (Use + (Ul ® Is)Uss(U s <8 h)) MO (u se - (U s ® h)U S£ (U s ® I s )) \i>)\£). (10) 
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(V) Measure the ancilla, thereby selecting the effective 
evolution on the joint Hs^Ti-e space. Such a conditional 
evolution will not be unitary in general. In particular, if 
the outcome is 0, then the conditional joint state may be 
seen as the result of the following generalized measure- 
ment operation: 

Pse(T) = -!- (a<^(T) Ps£ (0)A^(T)) , 

where p S£ {0) = \ip)(ip\ ® and 

A se = \ fee ± (Ul <8> h)Use{U s <8 Is)) , (H) 



po(T) = p(T)| o = Tr (4 + )t(T)A( + >(T)^(0)) , (12) 

is the corresponding conditional evolution probability at 
time T. From Eq. (|10f> . we notice that the net effect 
of the above (I)-(V) protocol is similar to a two-steps 
DD strategy [13, 123, where, however, the role of the 
unwanted coupling Hamiltonian is taken by the joint S£ 
propagator. 

From a control standpoint, one is naturally led to the 
following design issue: Given a joint propagator Use{T) 
for arbitrary T > 0, find unitary matrices Us and {/ fb on 
the system, such that the conditional decoupling condi- 
tions are satisfied, 



equivalent to 



a£(t) = uI®b ( £ -\t), 



(+)/ 



(13) 
(14) 



where -Bg (T) account for the effect of the dynamics on 
£ , and U fh conditional to the outcome |1)^ corrects the 
unintended component. Notice that if condition l|13(l can 
be satisfied with B £ unitary, the probability of obtain- 
ing 1 0)_4 is automatically one. In practice, although no 
constraint on the evolution time exists, ensuring perfect 
knowledge of the joint propagator Use{T) may be hard 
to achieve. We arc thus brought to determine conditions 
on Uss{T) that may still allow a solution of the above 
DD problem. 



A. Short time scale: Necessary and sufficient 
conditions 

Assume first that the relevant evolution time is suffi- 
ciently short (formally infinitesimal), so that a first-order 
Taylor expansion accurately represents the joint propa- 
gator Use, 



Use « Ise - iH S£ 5t + 0(St 2 



(15) 



The "decoupling-like" conditions we are seeking may be 
derived by substituting in (fHj)l . Due to the trace- 
lessness assumption on Si,i > 0, the latter equation is 



(Si + U s S. t U s ) ®B l= 



i>0 



The following results, which extends [23j, is relevant. 

Theorem 1 . Let X be an n- dimensional, traceless nor- 
mal matrix and sp(A) = {xi}™ =1 its spectrum. Then, the 
equation 



X + u f xu = , 



(16) 



admits a unitary solution U if and only if there exists 
a permutation of the elements of sp(X) such that the 
following mixing condition is satisfied: 



1. 



(17) 



//, in addition, X is Hermitian, and the spectrum el- 
ements are taken with descent ordering, then condition 
(II 711 implies a solution of the form 



U = VJV^ , 



(18) 



where V is a unitary matrix satisfying V^XV = X = 
diag(a;i, . . . , x n ), and 



J = 















Proof. First, notice that Eq. (|16fl has a unitary solution 
U if and only if 



X + ffJtU = 



(19) 



has a unitary solution U (in which case, U = V'UV). 
For the if implication, assume that (|17f) is satisfied for 
a permutation of the first m integers . . . , i m ). Then 
the unitary matrix U which reorders the basis vectors 
indexes according to that permutation, i.e., 

u \ipj) = \i>ij)> J = !)•••>'", 

satisfies (|19f) . For proving the opposite implication, as- 
sume that U is a solution of (EH). Then U^XU must be 
diagonal. Moreover, since a unitary change of basis does 
not change the spectrum, U'XU = diag(xi 1 , . . . , Xi m ), 
where . . . , i m ) is a permutation of the first m integers. 
Employing again (|19|) . we get X = — diag(a;, 15 . . . , Xi n ), 
hence condition (|17p holds. The specialization to the 
Hermitian case is straightforward. □ 

The above theorem fully characterizes the interactions 
able to be suppressed via a two-steps DD scheme: In fact, 
all the operators Si, i > need to satisfy the condition 
of the theorem, in the same basis, that is, employing 
the same unitary change of basis V. Moreover, the form 
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of the control actions to be employed is constructively 
provided in the proof. 

Assume that the relevant conditions are satisfied. 
Then we can choose Us such that 



A { /J(5t)=I s ®(I £ -iH s 6t), 
Thus, in this limit, 



i>0 



lim p(St) 
st^o 



0(St 2 



0, 



(20) 



(21) 



meaning that the faulty outcome 1 has zero probability 
to be obtained - consistent with the fact that B^(dt) 
is unitary. 

In principle, evolution over a finite time T may always 
be obtained by iterating the the basic steps (I)-(V) a 
sufficiently large number of times N, so that St = T/N — > 
0. The joint probability of obtaining zero as result of all 
the measurements up to T becomes then 



P(T) 



= lim 



N 



\m) 



lim 



N- 



-iH £ T 



i.T \ N 

wif>ir=i. 



as expected. In this form, our FDD protocol may be 
thought of as achieving an implementation of DD di- 



rectly based on the quantum Zeno effect |2j|. Formally, 
the protocol may also be regarded as a special case of 
the dynamical symmetrization scheme via repeated mea- 
surements originally proposed by Zanardi [13. While 
we have intentionally restricted the analysis to a two- 
dimensional ancilla (hence group-algebra in the language 
of this allows an explicit characterization (via The- 
orem 1) of the open-system Hamiltonians for which the 
protocol is effective. In addition, as it will be clear in the 
next Section, one of our main goals is to explore the use 
of classical feedback blocks to enhance the performance 
of open- loop DD. 

Before moving to consider finite-time feedback evo- 
lutions, we address the robustness of the FDD protocol 
against unintended Hamiltonian drifts on the ancilla sub- 
system, H a ^ 0. By starting from the same initial state 
for 5.4 as in step (I) above, in addition to the joint 
system-environment evolution Use, we now a ls° consider 
a unitary evolution U a on A. It suffices to describe its 
action on the basis states: 



|0) i-> cos0|O) + sin0e^|l) 
|1) 1 — * -sin0|O) +cos0e-**|l) 



(22) 



where 9, <f> are linear functions of t. Thus, the joint state 
after step (II) is 

W) = (cos0|O)+sin0e^|l})[/ ss |V>}|O 

+ (-sme\0} + cosee-^\l))Uss(Us\ip})\0- 
Applying C^ t and Ug as in steps (III)- (IV), respectively, 
yields: 



W) = |0) (cosOUse + sin0e^(£/J ® I £ )U S £ - sm9U S£ {U s ® h) + cosfcr*^ ® I £ )U S£ {Us ® I £ j) 
+ |1) (cos6U S£ - smee l4, (U s ® I £ )U S£ - sin9U S£ (U s ® h) - cos9(U ] s ® I £ )U S£ (U S ® h)) 

I 



Finally, if the joint probability for obtaining always |0) as 
a result of the measurements is computed, to first order 
in t one obtains 



P(T) 



lim 



0,...,0 



lim 



I--(H £ 
N 1 



, 0o , (Ul-Us) 



2i 



N 2 



where 0o = 0(0), #0 = #(0), respectively. Accordingly, 
if Us is both unitary and Hermitian (as it is the case 
e.g., for a CNOT gate), the same expression of the ideal 
case is recovered, and the above probability still tends 
to 1. Since, under these conditions, the unitary drift of 
the ancilla is irrelevant, this relax the requirement of a 
"frozen" ancilla for sufficiently fast protocols. 



B. Finite evolution time: Necessary and sufficient 
conditions 



We now provide necessary and sufficient conditions on 
the joint S£ propagator for employing the FDD protocol 
one-shot over a finite time interval. We begin with a 
definition. 

Definition 1 . An operator Xs £ acting on a bipartite 
Hilbert space TLse = T~ts ® T^E is locally diagonalizable 
(LD) if there exist operators {Aj\j and {Bj}j on TLs and 
He, respectively, such that 



Xs £ = Aj ® Bj , with 
3 

[A i ,A j ]=Q, Vi,j. 



(23) 
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The fact that an operator on a tensor product space 
can always be decomposed in the form Xgg = . Aj®Bj 
follows for example from the operator-Schmidt decompo- 
sition j2^| . Condition H23|) ensures that there exist a uni- 
tary change of basis, corresponding to an operator Us in 
TCs, that simultaneously diagonalizes all the Aj's, that 
is, Us brings Use to a block-diagonal form. 

Clearly, all factorizable normal operators on finite- 
dimensional Hilbert space are LD. However, the converse 
is not true. For Hamiltonians as in Eq. all evolutions 
generated under the condition that [Hs,Hss] = are 
LD. This includes, for instance, purely dephasing spin- 
bath or spin-boson models as in |29( . In the language 
of quantum error correction, LD open-system evolutions 
correspond to purely Abelian error algebras [30l | . 

Proposition 1 . Let dim(7Ys) = 2, and let the joint 
evolution with the environment £ at time T be described 
by the propagator Uss(T). There exists a one-bit feedback 
DD protocol that restores the initial state of the system 
for arbitrary T if and only if Use is LD. 

Proof. For convenience, assume now the ordering 8 (g) S 
in the tensor products. By hypothesis, Use satisfies |J23J|, 
thus we may work in a basis where all the Aj are diagonal 
by applying a suitable Us- In such a basis each block of 
U'se = UsUseUg becomes diagonal. Consider then these 
diagonal blocks, that is, the dg 2 x 2-matrices 









^2 ! 



I = 1, 



Each of them may be rewritten as: 









diA 



1>2,i 



,,d £ 2 . 



d\A 



Hence, the following decomposition applies: 

U'ss = B' 1 ®I 2 + B' 2 ®a z , 

where B[ , B' 2 are operators on Tig . It is now apparent 
from Eq. l(T3|l that, for example, Us = Uscr x U s can de- 
couple the system if the ancilla is found in the |0) state, 
Us being the change of basis that diagonalizes the Aj . If 
the outcome |1) is measured, the resulting evolution op- 
erator for the system and the bath in the diagonal basis 
is of the form: 

2B' 2 ®a Zl 

that is, factorized and unitarily correctable in the original 
basis for S by employing U lb = Uso- z U s . Necessity of the 
LD condition follows from the fact that each block must 
be proportional to the identity simultaneously and from 
invoking Theorem ^ □ 

For dimension, ds > 2, the algebraic conditions that 
need to be satisfied become more demanding: Essentially, 
the same kind of mixing structure found in Theorem ^ 
has to hold for Use- In fact, the above Proposition may 
be seen as a specialization of the following result: 



Theorem 2. There exists a one-bit feedback- DD proto- 
col that restores the initial state of the system for arbi- 
trary T if and only if Use is LD with each ds -dimensional 
block of the form: 



/d hi ••• \ 
d 2l : 



V 



'•• 

• • • d ds4 J 



ds 



L d s 



c t D, 



(24) 

where Cj is a coefficient depending on the block index, 
and D is unitary and satisfies the spectral conditions of 
Theorem 1. 

Proof. Since Use is LD and Eq. holds, the fact 
that D satisfies (I17f) is a sufficient and necessary condi- 
tion for the existence of an effective Us by use of Theorem 
n Here, however, we deal with complex eigenvalues, so 
there is no descent ordering and we cannot give the form 
of decoupling operation Us in a natural way. It has now 
to be proven that there exists a unitary operation work- 
ing for the feedback correction step if the faulty outcome 
has been measured. From Eq. 114(1. this is true if and 
only if D is proportional to a unitary operation, i.e. all 
its eigenvalues have the same absolute value. Including 
the coefficient in c, yields the desired result. □ 

The FDD strategy has a potentially important advan- 
tage with respect to the standard open-loop averaging 
schemes, provided that the appropriate algebraic con- 
ditions are fulfilled: Its effectiveness does not directly 
depend on the elapsed interaction time T. In practice, 
of course, several characteristic time scales, depending on 
both the open-system and control operations, will still be 
relevant. In particular, we assumed to be able to enact 
ideal, instantaneous unitary preparation, measurement, 
and resetting steps throughout. In practical terms, this 
will require the capability of implement fast conditional 
gates and measurements with respect to the typical cor- 
relation time scale of the underlying open-system dynam- 
ics. The other strong assumption we made is to be able 
to access an effectively decoherence-free ancilla |3lj . 



IV. FEEDBACK-CONCATENATED 
DYNAMICAL DECOUPLING 

In the previous Section, we developed necessary and 
sufficient linear-algebraic conditions for the unitary, joint 
evolution operator for the target system plus environ- 
ment Use to be fully correctable by a feedback strategy 
involving a combination of fast conditional operations 
and projective measurements. However, making sure 
that such conditions are verified, and explicitly construct- 
ing the required control operations, demands in general 
an accurate knowledge of both the form and the dura- 
tion of the interaction - restricting the usefulness of the 
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procedure in realistic scenarios. If such a detailed infor- 
mation is not available, our next goal is to show how one 
can still make a feedback strategy effective by combining 
the feedback action with open-loop averaging schemes - 
in particular, by appropriately modifying the evolution 
between subsequent measurements. 

In the rest of our present discussion, we will develop a 
general solution in the simple yet paradigmatic case of a 
single qubit (d$ = 2) as the system of interest. Possible 
extensions to higher-dimensional systems will be men- 
tioned in the conclusions. 



A. Problem setting 

For the sake of simplicity, we will describe in detail the 
basic strategy for the Hamiltonian suppression case, that 
is, we set Bi = for all i in Eq. QJ. Notice that this limit 
could still account for non-unitary effects from a scmiclas- 
sical environment if randomly fluctuating parameters are 
allowed in the Hamiltonian - prominent physical exam- 
ples being decoherence from random telegraph noise |32l 
and from a dipolarly coupled nuclear spin reservoir |33j| . 
The proof for the general open quantum system case re- 
lies on Proposition^ The control problem for which we 
will provide a solution may then be stated as follows: 

Problem 1 . Consider a single qubit in the state p(0) at 
t = 0, undergoing an uncertain Hamiltonian evolution. 
Assume we have an error set V of possible Hamiltonians 
{H}, and an initial estimate H G V for H. The task is 
to ensure that p(T) = p(0) at a given finite time T , for 
any possible H G V. 

This problem, which is equivalent to the implementa- 
tion of a no-op unitary gate between times and T, qual- 
ifies as a robust control problem. That is, the output state 
must be insensitive to errors in H . The most general error 
Hamiltonian may be described as an element of i ■ su(2), 
that is, AH = H — H = e- cr = e x a x + e y <J y + e 2 cr z . The 
initial estimate H is given by some a priori knowledge, 
or guess, on the system's behavior. It is often the case 
that approximated models are available for the dynamics, 
and if such additional information is available, a natural 
question is whether and how such information may be 
efficiently exploited to improve over strategies which do 
not. In our case, a robust solution which invokes no an- 
cilla and feedback loop, and does not require in principle 
prior information, is provided by the open-loop, maximal 
DD scheme described in Section II [fT]]- In fact, for both 
purely open- loop and feedback-based strategies, the role 
of available prior information becomes critical as soon as 
ideal control conditions are relaxed. This will be quanti- 
tatively addressed in Section Ivl 



B. Concatenated control protocols 

Let as above H = H + AH represent the target qubit 
Hamiltonian in terms of an estimated plus error com- 
ponent. Consider a choice of axis and units such that 
H = a z . A first, natural way to merge open- and closed- 
loop DD is to employ a selective DD to remove the main 
component of the Hamiltonian, while relying on feedback 
to counteract the residual ones. We will refer to the re- 
sulting protocol as Feedback-Enhanced Decoupling (FED). 
From Section^ a singlc-qubit selective DD scheme based 
on the control propagator sequence {I, cr x } projects onto 
the a x direction in the ideal limit of T c — > 0. Thus, to 
first order, only an error component e x o~ x will be left on 
the estimate. The idea is to apply an "outer" block of 
FDD, so that a suitable choice of conditional operator, 
e.g., will correct for the above residual component. 

In analogy to steps (I)-(V) in Section lTTTl the FED pro- 
tocol proceeds as follows: 

(I') Prepare the ancilla in the initial state \<p) = 
(|0) + \l))/y/2 and entangle it with the qubit system by 
applying the operation. 

(IF) Evolve under a cycle of selective DD {I,a x }, by 
obtaining at T = T c a net propagator of the form: 

Ur(T c ) = e-**" , 

where <f> = e x T c . 

(III')-(IV') Undo the preparation steps applying the 
conditional gate and a Hadamard U^, respectively. 

(V') Measure the ancilla in the {|0), |1)} basis, obtain- 
ing outcome with probability po(T c ) = cos 2 (4>), and 1 
with probability pi(T c ) = sin 2 (0). 

Ug es {T c ) satisfies the conditions of Proposition ^ In 
particular, it is easy to see that Us — o~z is a suitable 
operation for counteracting such a residual evolution, ir- 
respective of the value of <p. Hence, as discussed previ- 
ously, the classical information extracted from the ancilla 
measurement restores the system's state in the "0" case, 
while in the other case it indicates univocally the cor- 
rection needed for taking the system back to the initial 
condition. In fact, applying a fast resetting operation 
o~x ® a x when the measurement outcome is one, we re- 
store the initial condition of both iS and A together. Re- 
markably, this happens no matter how large is the error 
we are making in the prior, thereby producing a maxi- 
mal DD strategy in the ideal, fast control approximation. 
This also means that the strategy is effective for every 
Hamiltonian in V, satisfying the robustness requirement. 

Note that at this level the role of the initial estima- 
tion does not emerge clearly, suggesting another possible 
concatenation scheme in place of the above FED strategy. 
The basic observation is that nothing prevents us, in prin- 
ciple, from swapping the roles of the feedback loop and 
the selective DD strategy. Indeed, we may use a {/, a z } 
selective DD sequence to obtain a net evolution of the 
form 

U/ es (T c )=e~^'^ , 
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TABLE I: Overview of the relevant control protocols for a single qubit, in the simple case of Hamiltonian suppression. We 
assume H — H + AH, where H is the Hamiltonian estimate. NC stands for "Not Corrected" . 



Protocol 


Acronym 


H corrected by: 


AH corrected by: 


Selective Decoupling 


SelDD 


Decoupling 


NC 


Maximal Decoupling 


MaxDD 


Decoupling 


Decoupling 


Feedback-enacted DD 


FDD 


Feedback 


NC 


Feedback-Enhanced Decoupling 


FED 


Decoupling 


Feedback 


Decoupling-Enhanced Feedback 


DEF 


Feedback 


Decoupling 



where now <f)' = (1 + £z)T c , and employ operations 
along with a a£ (g> a z correction for the feedback loop. 
This choice leads to another effective solution for our 
problem, in which we use the selective DD to remove the 
errors, and the feedback loop to correct H. We will refer 
to this second protocol as Decoupling-Enhanced Feedback 



(DEF). We will compare the effect of non-ideal condition 
for both strategies with the aid of numerical simulations 
in Section[V] A summary of the relevant strategies, along 
with the respective roles of the feedback and open-loop 
components, is given in Table [I] 



C. Estimation and tuning 

The FED procedure described above offers a robust so- 
lution to Problem ^ Nevertheless, it also provides us 
with additional information we have not used yet: The 
record of subsequent measurement outcomes. In fact, 
the feedback correction we apply needs only the last 
outcome, and no extra memory. However, by record- 
ing individual outcomes, we have (under suitable sta- 
tionarity assumptions) a way for estimating (e^l, i.e., 
\e x \ = arcsin(VpU%))/^c- 

On the other hand, if a DEF protocol is implemented 
instead, this allows us information on \e z \ to be acquired 
- and, with suitable rearrangements, on \e y \ as well. We 
will illustrate in the next Section how better informa- 
tion about the underlying Hamiltonian is important to 
optimize control performance when ideal conditions are 
relaxed. 

Thus, the acquired information may be exploited to 
tune the overall control strategy on a better estimate a. 
Ultimately, this may result into iteratively converging to 
the correct Hamiltonian. In turn, this will decrease the 
need for control actions in the feedback correction step, 
as well as the impact of higher-order errors induced in 
the concatenated FED/DEF protocols for finite T c . 

Attention must be payed to the fact that we are only 
gathering information about the absolute value of e x , z , 
and to the fact that we need to estimate the component 
e y in a second stage, by replacing a z with a y operations in 
the relevant DD protocol. In what follows, for simplicity 
we will not consider errors in the y direction. The main 
steps of a strategy for Hamiltonian estimation may then 
be sketched as follows: 

(i) Run the FED protocol M times. Consider a z as 
estimate and record in R the sum of the results of the 



measurement steps. Compute pf = R/M and the esti- 
mation for |e x | = arcsin(y/pf)/T c . 

(ii) Switch to the DEF protocol, and run it M times. 
Compute p\ = R/M and the estimation for e z = 
(arcsin(y / pf) - 1)/T C . 

(iii) Update the estimate, a z ew = a z +i] z e z a z ±r] x e x a x , 
where < r\ z ,r\ x < 1 are two weighting parameters to 
tunc, depending on the significance of the estimate for 
finite M. In the first iteration, the sign of r\ x can be 
chosen arbitrarily. 

(iv) Iterate steps (i)-(m), changing the sign of i] x if 
the frequency of corrections pf is increasing - indicating 
a worse estimate. This shall complement on average the 
sign information missing on e x . 

The above procedure will approach the best approx- 
imation of the coupling Hamiltonian in the x-z plane, 
compatibly with the estimation errors due to a finite 
number of trials (in any non ideal setting, i.e. At > 0, M 
is finite). At a subsequent stage, it is certainly possible to 
switch to a FED strategy employing o~ y transformations, 
and refine the estimation in the x-y plane also. 

V. NUMERICAL RESULTS 

If, theoretically, both maximal and feedback-enhanced 
DD attain the desired dynamical averaging, from a prac- 
tical standpoint it is critical to assess their performance 
as the ideal conditions are relaxed. In particular, we will 
focus on comparing the response of different single-qubit 
control strategies when the distance At between two ideal 
pulses remains finite, and we will elucidate the role of the 
prior information on the Hamiltonian. 

As a state-independent performance indicator for con- 
trol performance, we invoke entanglement fidelity T |34j| . 
Borrowing from standard information-theoretic results 
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for singlc-qubit quantum channels, the latter will be com- 
puted indirectl y fro m the average fidelity or from gate fi- 
delity, see e.g. 

jus HIS Wewi11 

use, in particular, 

the useful result that determines the average fidelity for 
a generic quantum operation T as: 

P{T) = l + T2 ^ trace fo-Tfo)). (25) 

Then, T may be explicitly computed as 

d 

where d = ds = 2 in our case. Additionally, if T is 
unitary, T(aj) = UajU^ , T(T) simply reduces to 

T(T) = ^ |trace(t/(T))| 2 . (26) 

In order to ease the comparison, let us introduce some 
additional compact notation for the various strategies un- 
der investigation. We will consider two selective strate- 
gies, i.e. strategics that are meant to remove only the 
estimated Hamiltonian component: The selective, sym- 
metrized two-pulses Carr-Purcell sequence described in 
Section|n] (SelDD), and the FDD strategy introduced and 
discussed in Section IIIII Next, we will compare some 
maximal DD schemes, i.e. strategics that in the ideal 
limit At — > remove any Hamiltonian component of the 
evolution: The maximal, four-pulses open-loop DD se- 
quence sketched in Section [H] (MaxDD), and the two con- 
catenated strategies presented in Section HVI the FED and 
DEF protocols. 

Before discussing in detail the numerical results, it is 
worth anticipating that a clear asymmetry emerges in the 
observed behavior with respect to errors on the estimated 
Hamiltonian. For example, the FED and DEF performance 
changes if the faulty Hamiltonian employed in the simula- 
tions is a z +ea x or a z + ea y , respectively (a z representing 
the estimated Hamiltonian, as in the previous Section). 
In a similar fashion, different, but still ideally equivalent 
pulse sequences for MaxDD, result in non uniform perfor- 
mance with respect to a given Hamiltonian. 

Higher-order errors emerging from the Magnus series 
are responsible for such asymmetry. If H is the Hamil- 
tonian wc wish to suppress with a sequence {Gi}™=i °f 
unitary control actions, define Hi = G[HGi. Then, re- 
calling Eq. ©, the first order correction to the average 
Hamiltonian in © rewrites as 

Thus, it is straightforward to sec that, for instance, for a 
Hamiltonian H = a z + sHi, with trace(er z i/i) = 0, the 
control sequence {/, a x , a z , a y } in Eq. 127|) gives non-zero 
terms only to second order in e, 

B{ X) = ^ £ 2 [axHi<7xjHil (2g) 



whereas the sequence {I,a y ,a x ,a z }, gives non-trivial 
terms already to first-order in e: 

o - Aj.2 

s(1) = fl^ e[(JyHi(Jy + Huax]+o{£ 2y (2Q) 
-t c 

Note that in both cases the corrections are zero for e = 0, 
showing how, could we know the exact Hamiltonian and 
choose the pulses freely, we could optimize control per- 
formance. Thus, even simple DD strategies may signifi- 
cantly benefit, in general, from a better knowledge of the 
Hamiltonian. (In fact, accurate knowledge of the under- 
lying Hamiltonian is one of the key elements for high-level 
DD design, as exemplified by a large variety of multipulse 
DD sequences in NMR jH). 

On the other hand, for a selective DD scheme based 
e.g. on {/, a x } and the same Hamiltonian H = a z + eH\ 
as above, to first order in e wc obtain 

H (1) = — ^—[crxHicr x + Hi, a z ] , (30) 

which is proportional to a y . This shows how, in order to 
retain only the z-component of the Hamiltonian with a 
{/, er 2 } scheme, the corrections due to a non-zero At will 
be along y. This means that a feedback strategy employ- 
ing would correct them along with the z-component 
of the Hamiltonian, while the use of would not. Such 

y 

a perturbative analysis is consistent with the results of 
the simulation, and accounts for the observed asymme- 
tries of both the maximal DD and feedback based strate- 
gies. These observations offer useful insights on the best 
choice of control actions once the main component of the 
Hamiltonian is known. 

We now proceed to illustrate and discuss some repre- 
sentative numerical results, for a generic qubit Hamilto- 
nian of the form H = £l z a z + e x a x + s y a y . First, we 
compare the SelDD with the FDD protocol in the presence 
of a finite intcrpulsc separation At (which is the control 
non-ideality we focus on). From Eq. 112611 . the free evolu- 
tion corresponds to entanglement fidelity 

F{T) = cos(OT) 2 , Q w Q z . 

A typical behavior is depicted in Figure ^ The changes 
in the oscillation frequency, Q e x , i— > £ y , associ- 
ated with each SelDD protocols, are clearly seen. Not 
only does FDD display better performance for compara- 
ble, fixed At, but it also turns out to be less sensitive to 
errors in the estimate of the Hamiltonian. In fact, the 
FDD protocol exhibits better performance as the norm of 
the estimation error grows. 

In comparing maximal DD strategies, we made sure 
to employ in each case the appropriate choice of control 
pulses, as following from the analysis of the higher-order 
corrections presented above. Different control scenarios 
were investigated. The main features emerging from our 
study may be summarized as follows. 

As the parameter At grows within the parameter range 
considered, all protocols move away from the ideal behav- 
ior in a regular fashion: Overall, the best performance is 
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FIG. 1: Comparison between selective DD protocols: 

(color online) FDD (brown, circles) and SelDD, with a x (blue, 
squares) or a y (green, triangles) pulses. Symbols are displayed 
every four decoupling cycles to improve readability. Brown 
and blue curves correspond to oscillations of the form T(T) = 
cos{e x , y T) 2 , respectively - as following from Eq. (1261 for the 
selectively decoupled evolution. Actual Hamiltonian: H — 
£l z o z + s x a x + Ey&y, Ex = fl z /20, £ y = f2,/10. The estimated 
Hamiltonian is proportional to a z . flT = 30, Q z At = 0.04. 



FIG. 2: Comparison between maximal DD proto- 
cols: (color online) MaxDD (blue, squares), symmetrized FED 
(brown, circles), non-symmetrized FED (green, triangles), DEF 
(red, diamonds). Symbols are displayed every two decou- 
pling cycles to improve readability. Actual Hamiltonian: 
H = Q z a z + e x a x + e y a y , e x = e y = fl z /10. The estimated 
Hamiltonian is proportional to a z . Q Z T = 120, Q z At = 0.32. 



VI. DISCUSSION AND CONCLUSION 



attained by the FED strategy, followed by MaxDD. The DEF 
strategy, where we use the feedback loop to correct the 
main, estimated Hamiltonian, shows a sensible reduction 
of the performance for longer At, making MaxDD more 
effective in the long time regime. Independent tests con- 
firmed that this implementation is more sensitive to the 
At parameter increasing. For sufficiently small At, DEF 
can outperform the maxDD, however, compared to FED, it 
looses its advantage soon. 

For a fixed value of both At and T, we also investigated 
the performance of the above protocols when the norm of 
the estimation error on the Hamiltonian is increased. The 
same comparative results are found, with fairly uniform 
behavior, and no particular crossings. 

Based on the above analysis, the most robust and 
information-efficient strategy turns out to be the sym- 
metrized FED protocol, that is, the feedback-corrected 
version of a symmetric selective DD. It is worth to notice 
that the symmetrization on nested SelDD is crucial in 
determining such advantage: The non-symmetrized ver- 
sion does not succeed, in general, at outperforming the 
open- loop strategy. We summarize in Figure [21 the typ- 
ical behavior of the protocols under considerations for a 
fixed At value. We show also the simulation results for 
a FED strategy with a non-symmetrized nested SelDD, 
clearly supporting the above claim on the importance of 
symmetrization. 



We have designed and characterized different quan- 
tum procedures for exploiting single-bit classical feedback 
as a tool for quantum dynamical averaging - applica- 
ble alone or in conjunction with traditional (determin- 
istic) open-loop decoupling methods. Beside providing 
general (linear-algebraic) necessary and sufficient con- 
ditions for a two-step averaging processes to be imple- 
mentable in principle via single-bit feedback on an arbi- 
trary finite-dimensional open quantum system, our anal- 
ysis points to novel possibilities for successfully merging 
open- and closed-loop techniques toward achieving robust 
quantum dynamical control. In particular, the concate- 
nated feedback-enhanced decoupling protocol proposed 
and discussed for a two-dimensional target system in Sec- 
tions IIVIVI is found to exhibit better performance than 
its open-loop counterpart, the best possible maximal DD 
scheme with respect to the estimated interaction. For 
higher dimensional systems, the integration of discrete- 
time feedback with open-loop decoupling schemes ap- 
pears more delicate, as the constraints imposed by The- 
orem |2 become more demanding. While it is still con- 
ceivable in principle to employ feedback to enhance the 
performance of selective open-loop strategies, especially 
whenever uncertainty is limited or it affects the free evo- 
lution in some structured fashion, further investigation is 
needed to address the relevant issues in detail. 

From an information-theoretic perspective, our results 
offer some suggestive insights on the role played by clas- 
sical information in various quantum control schemes. 
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While different schemes become essentially equivalent 
and insensitive to system's parameters in the ideal limit 
of arbitrarily fast control, prior information on the un- 
derlying Hamiltonian plays a central role in the optimiza- 
tion of different protocols as soon as parameter uncer- 
tainties become important and non-ideal control settings 
are considered - as it is unavoidably the case in prac- 
tice. Thus, the role of information in control design and 
synthesis cannot be underestimated. Interestingly, even 
if purely open-loop maximal decoupling schemes may be 
optimized by using the initial estimation, they do not al- 
low any improvement in the estimate itself. Thus, the full 
power of a feedback strategy might is expected to appear 
in adaptive control settings, like the simple Hamiltonian 
estimation setting presented in Section HV CI 

From yet another standpoint, the proposed feedback- 
concatenated protocols further support the relevance of 
hybrid strategies for robust quantum dynamical control 
|39| . In particular, our results take an additional step to- 
ward establishing the potential of dynamical decoupling 
methods for integration with existing passive and active 
quantum control techniques - following the demonstra- 
tion of concatenated quantum error correction and de- 
coupling codes (ZJj, the development [2f| 0, ^3] an d 
implementation of dccohcrence-free encoded dynamical 
decoupling schemes |37j . and the recent development of 
fault-tolerant, concatenated dynamical decoupling proto- 
cols 01 • Beside, as mentioned earlier, generalizations of 
the proposed framework to higher-dimensional systems 
and ancillas, additional extensions to include non-ideal 
measurement and realistic control capabilities are cer- 
tainly necessary for a more complete picture and possi- 
ble contact with experiment to be established. At the 
formal level, we expect that linear-algebraic tools may 
continue to prove useful in that respect, along with gen- 
eral methods from quantum fault-tolerance theory. In 
practice, although meeting all the relevant control re- 
quirements appears demanding by present capabilities, 



the kind of operations needed for the proposed hybrid 
feedback-decoupling protocols have already been imple- 
mented separately in different device settings Q ( see 
also discussion in |14|L and experimental progress is 
steady. In this respect, it might be especially intriguing 
and timely to explore possible proof-of-principle appli- 
cations in scalable devices based on trapped ions, where 
quantum error correction has been recently demonstrated 
[45j and different isotopic species could be used i n p rin- 
ciple for system and ancilla degrees of freedom |4(j, or 
semiconductor quantum dots, where suggestive propos- 
als for controlling a slowly fluctuating mesoscopic spin 
environment via a suitable probe spin already exist [4 71 ]. 

Lastly, the analysis of higher order corrections, as car- 
ried out in Section may be seen as a simple instance 
of the sensitivity minimization problem, which is one of 
the fundamental problems of robust control theory. In 
this perspective, it is our hope that the present work 
will stimulate further exploration and input from classi- 
cal control theory, as well as the quantum control and 
system engineering communities. 
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